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An ideal material for on-board hydrogen stor- 
age must release hydrogen at practical tempera- 
ture and pressure and also regenerate efficiently 
under similarly gentle conditions. Therefore, 
thermodynamically, the hydride material must 
lie within a narrow range near the hydrogena- 
tion/dehydrogenation phase boundary. Materi- 
als involving only conventional bonding mecha- 
nisms are unlikely to meet these requirements. 
In contrast, materials containing certain frus- 
trated bonding are designed to be on the verge 
of frustration-induced phase transition, and they 
may be better suited for hydrogen storage. Here 
we propose a novel layered solid boron hydride 
and show its potential for hydrogen storage. The 
absence of soft phonon modes confirms the dy- 
namical stability of the structure. Charging the 
structure significantly softens hydrogen-related 
phonon modes. Boron-related phonons, in con- 
trast, are either hardened or not significantly af- 
fected by electron doping. These results suggest 
that electrochemical charging may facilitate hy- 
drogen release while the underlying boron net- 
work remains intact for subsequent rehydrogena- 
tion. 

Hydrogen fuel cells are attractive for their very high 
efficiency. Beyond the technology of the cells themselves, 
these systems present two major challenges to current 
technologies: economic large-scale hydrogen production 
and practical small-scale on-board hydrogen storag 
The storage problem in particular demands fundamen- 
tally new transformative science and engineering, since 
the standard modes of chemical and physical interac- 
tion known for hydrogen appear to be grossly insufficient 
to meet the needs. All current hydrogen storage tech- 
nologies (e.g. cryogenic liquid, high-pressure gas cells, 
low-temperature adsorbates, metal hydrides, and chem- 
ical storage) suffer from one or more major deficiencies 
in cost, thermodynamic efficiency, volumetric capacity, 
kinetics, gravimetric density, long-term storage and re- 
versibility^. These deficiencies are so strong that a clear 
engineering path to a practical storage system has not 
yet been envisioned. For example, transition metal hy- 
drides usually have desired hydrogen release kinetics but 
fall short of the target storage capacity. Physisorption 
storage schemes, on the other hand, suffer from thermo- 
dynamic issues of storage capacity. Complex hydrides 



formed from the restricted class of light metals (e.g. 
MA1HJ3 and MBHj^, where M denotes light metal such 
as Na or Li) and amide materials^ suffer from issues of 
kinetics and variable desorption energies for successive 
hydrogens release. Systems based on amino-boranes (e.g. 
NHsBHjE^ and related compounds) have a fine balance 
of chemical energies that enable hydrogen release under 
relatively gentle conditions. However, successive hydro- 
gens are progressively harder to remove, and the final 
dehydrogenated compound is a refractory ceramic that 
is not amenable to reversible rehydrogenation. 

It is interesting that all three major classes of hydro- 
gen storage materials investigated so far (i.e., low tem- 
perature adsorbates, metal hydrides, and chemical stor- 
age) overlaps with regards to only boron. Examples 
include physical adsorption systems with boron-doped 
carbonsSHISl, chemical systems using ammonia borane 
derivatives and light metal borohydrides. Boron is dis- 
tinguished from the remainder of the periodic table by an 
electronic frustration. As a result, boron-based materials 
are well-known for their rich and highly unconventional 
chemistry. Having only three valence electrons, boron 
cannot form a sufficient number of traditional covalent 
bonds to reach a close-shell configuration. Being a small, 
tight first-row p-block element, it is insufficiently elec- 
tropositive to participate strongly in traditional metallic 
bonding, yet also insufficiently electronegative to be dom- 
inated by highly stable ionic configurations. These prop- 
erties force boron to assume unusual multi-center bond- 
ing configurations in many of its compounds. One of the 
most famous examples, and a distinct feature of electron 
deficiency, is the three-center two-electron (3c2e) found 
in many boranes. Chemically frustrated boron, when 
properly deployed in a supporting chemical framework, 
provides an intriguing opportunity to attain an interme- 
diate bonding with hydrogen. In this paper, we exploit 
the rich chemistry of boron and hydrogen to propose a 
novel solid boron hydride structure which contains inter- 
mediate strength 3c2e-like multi-center B— H bonds and 
investigate its potential for on-board hydrogen storage. 



Model 

Although the rich boron chemistry explored to date 
has been dominated by three dimensional cage struc- 
ture^UM^, computational and experimental studies^H^ 
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have begun to uncover a latent planar aspect in boron- 
rich materials. In the past decades, ab-initio results 
of different planar and quasi-planar structures of boron 
clusters have been proposecP^^. Though these clus- 
ters were studied experimentally, understanding of their 
structural and electronic properties has been a challenge. 
Most rece ntly, photoelectron spectroscopy (PES) tech- 
niqu o 20 ! 21 ! is used to unravel better structural informa- 
tion of boron clusters. Using a combination of PES tech- 
nique and ab-initio calculations, a stable planar struct ure 
of boron clusters containing ten to fifteen atonis^l2!J are 
shown to exist. Furthermore, a new class of planar boron 
sheet s are predicted to exist through a self-doping mech- 
amglr j22i23] rp^g g^ajjjifty f these sheets are explained 
with a bonding mechanism which arises from the compe- 
tition of a two- and three-center bonding. 

Pure hexagonal boron layer, although not stable alone 
due to electron deficiency, can be stabilized by interca- 
lating electron donors such as Mg and Al to form MgB2 
and A1B 2 , (and also TaB 2 , HfB 2 , ZrB 2 , and TiB 2 ). This 
hexagonal boron sheet has a much greater structural 
flexibility than does graphite: for example, the in-plane 
lattice constant expands from 3.01 A in A1B 2 to 3.17 A 
in ZrB 2 . Very recently, the possibility of using planar 
boron sheets with possible alkali metals as a template 
for hydrogen storage has been investigated^. In this 
work, we propose that the electron donor can also be 
hydrogen atoms. In fact, nature finds a way to sta- 
bilize the electron-deficient molecule diborane through 
the formation of 3c2e bonds between boron and hydro- 
gen. It is plausible that the electron-deficient hexagonal 
boron network can be stabilized in part or whole through 
the formation of 3c2e bonds with hydrogen, as shown in 
Fig. [T] which we denote as "B 2 H 2 " . The boron atoms in 
B 2 H 2 form borane-like 3c2e bonds with bridging hydro- 
gen atoms and connect to their third in-plane neighbor 
via a conventional sp 2 hybridization (shown in Fig. [lj). 
Note that the hydrogen atoms are bound to the hexago- 
nal boron plane in pairs, one above, and the other below 
the plane. These weakly bonded hydrogens are crucial to 
stabilize the boron network and are desirable for the re- 
lease of hydrogen at practical conditions. The extended 
boron framework in this solid boron-hydride structure 
is quite distinct from a collection of individual diborane 
molecules. Preservation of the backbone boron struc- 
ture over cycling (by means of kinetic barriers against 
framework decomposition) would be extremely advanta- 
geous for regeneration as compared to traditional boron 
hydrides, i.e., boranes. 



Methods: Computational details 

The density functional theory^ based first-principles 
calculations in the present work are performed using 
the QUANTUM ESPRESSC^ package. The general- 
ized gradient approximation (GGA) of Perdew, Burke, 
and Ernzerhof 271 is used for the exchange and correlation 




FIG. 1: Structure of single- layer B2H2: a, The relaxed 
structure, b, The corresponding bond angles and 
bond lengths, c, The structure of diborane. Boron is 
blue and hydrogen is orange. 

terms. Ultrasoft pseudopotentials^ are used to describe 
the electron-ion interactions. The plane- wave energy cut- 
off is set at 50 Ry for the wave functions expansion and 
500 Ry for the charge density. All structures are relaxed 
until the forces on each atom are less than 10 -5 Ry/a.u. 
and the stress less than 1 GPa. In this work, we focus 
on the properties of single-layer B 2 B 2 . Inter-layer bind- 
ing will be investigated in the near future using a re- 
cently developed method for efficient calculations of the 
van der Waals energy^. The charge density for a single- 
layer B 2 H 2 is sampled with a fc-point set of 32x32x8 
in the Brillouin zone. The phonon dispersion is calcu- 
lated using the density functional perturbation theory^. 
The dynamical matrices are calculated on a 6x6x4 q- 
point. Next, the force constant matrices in real space 
are obtained by Fourier-transforming the dynamical ma- 
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trices. Finally, the dynamical matrix at any g-point is 
obtained by inverse Fourier-transformation of the force 
constant matrices. In the following sections, we discuss 
the electronic, structural, and dynamical properties for 
the proposed solid boron hydride B 2 H 2 . 




FIG. 2: Band structure of a single-layer of B2H2. 
The Fermi level is at eV and shown with red dotted 
line. 
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FIG. 3: Projected density of states (PDOS) and total 
density of states (DOStot) in a single-layer B2H2. a, 
PDOS for B 2s and H Is, and b, PDOS for B-2p x , 
B-2p„ and B-2p z . The Fermi level is at OcV. PDOS 
units are states/eV/cell. 



Structural and electronical properties of single-layer 

B2H2 



The proposed B 2 H 2 layered structure has a D 2 h sym- 
metry which is the same as that of diborane. The cal- 
culated bond angles and bond lengths of B 2 H 2 are close 
to those in diborane. The boron— boron (B— B) distance 
of the 3c2e bond is 1.82 A, to be compared with 1.76 A 
in diborane. The sp 2 B— B bond length (1.76 A) is sub- 
stantially shorter than the 3c2e B— B bond. For compari- 
son, the theoretical B— B bond length in MgB 2 calculated 
within the GGA is 1.77 A. The hydrogen-boron (H-B) 
bridge bond length is 1.33 A. This is essentially the same 
as that in diborane and is substantially longer than the 
sp 2 H-B bond (1.19 A). The B-H-B bridge bond an- 
gle is 86°, to be compared with 83° in diborane. These 
hydrogen bridge bonds are crucial to stabilize the boron 
network. 

The band structure and projected density of states of 
the single-layer B 2 H 2 are shown in Fig. [2] and Fig. [3] re- 
spectively. It is interesting to compare planar B 2 H 2 to 
MgB 2 . Like MgB 2 , the system is metallic. However, 
unlike MgB 2 , the degeneracy of the p — a bands at the 
T-point in MgB 2 is now removed due to lowering of the 
symmetry from D$h to D 2 h- As a result, only one of the 
p — d bands crosses the Fermi level and the hole pocket 
near the T-point comes from the upper p — a band. The 
degeneracy of the p — ir bonding and antibonding states 
near the Fermi level is also removed. Similar to that in 
MgB 2 , the Fermi surface of B 2 H 2 also consists of both a 
p — (7 derived hole pocket and a p — 7r derived electron 
pockets. The projected density of states plot also shows 
significant involvement of the hydrogen s states in the 
bonding. 



Phonon dispersion and dynamical stability 

To investigate the dynamical stability of the proposed 
structure, we carry out phonon calculations using density 
functional perturbation theory^. The calculated phonon 
dispersion (shown in Fig.[4| reveals no soft phonon modes 
(i.e., imaginary phonon energies), thus verifying the dy- 
namical stability of the B 2 H 2 structure. The nine optical 
phonon modes can be classified using group theory anal- 
ysis: 



4[H-H(z)] 4 

B 2g [R-U(y)} 

B lu [R-B(z)] 



A 2 g [B-B(x)]+B lg [B-B(y)} 



Bl g [R-R(x)} 
B 2u [R-B{x)} 



£ 3 2 ff [B-B(z)]4 
B 3u [R-B(y)]. 



The two A g modes will mix so that the polarization 
vectors of the eigenmodes with the A g symmetry will in- 
volve both hydrogen and boron motions. We denote the 
mode with primarily hydrogen motion as A*. Similarly, 
the two Bz g modes will mix and we denote the one with 
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FIG. 4: Phonon dispersion: a, The calculated phonon 
dispersion for single-layer B2H2. The labels at F are 
the nine optical phonon modes, b, The corresponding 
phonon polarization vectors of the optical phonons at 
r . Boron is blue and hydrogen is orange. 



primarily hydrogen motion as B 3g . The phonon polar- 
ization vectors for the 9 optical phonon modes at V are 
also shown in Fig. |4|b) . 

The highest energy mode (A*) is a symmetric stretch- 
ing of H atoms perpendicular to the boron plane (i.e. 
along the z-direction), whereas the B\ u mode involves 
H— B vibrations in the z-direction. The next two phonon 
modes {Bt, g and £?2u) ar e H— H and H— B vibrations 
along the direction of the bridging bonds (x-direction). 
The B\g and A g modes are similar to the E g modes in the 
AIB2 structure and involve boron motions in the boron 
plane (x-y plane). Note that the energies of these two 
modes are very close to the energy of the E g phonon in 
A1B 2 but substantially higher than that in MgB 2 . The 
B2 g and B$ u modes involve H— H and H— B vibrations in 
the y-direction whereas the low energy B 3g mode consists 
primarily of boron vibration in the y-direction admixed 
with hydrogen motion in the x-direction. 

One key aspect of the present work is to assess the fea- 
sibility of the B2H2 structure for hydrogen storage. This 



requires not only a stable backbone structure that binds 
hydrogen with intermediate strength but also a structure 
that releases hydrogen under practical conditions. In the 
next section, we investigate the doping effects on the dy- 
namical properties of B2H2 and discuss their implication 
for dehydrogenation. 



Charge doping effects 

Utilizing B2H2 for hydrogen storage requires not only 
a stable boron framework but also a practical hydro- 
gen release mechanism. As discussed earlier, a planar 
hexagonal boron layer is not stable alone due to elec- 
tron deficiency. If additional electrons are introduced to 
the B 2 H 2 system (for example, by introducing Li), the 
boron network will be less electron deficient. As a result, 
the hydrogen bridge bonds should weaken with increas- 
ing charge doping. The hexagonal boron backbone struc- 
ture, on the other hand, should remain stable under elec- 
tron doping. Therefore, one may be able to control the 
strength of the hydrogen bridge bonds through electro- 
chemical charge doping. This, when combined with ther- 
mal activation, may facilitate dehydrogenation and/or 
rehydrogenation. 

Here we investigate this possibility by studying the ef- 
fects of doping on the zone center phonon energies. Our 
results reveal that some of the phonon modes (especially 
those involving vibrations of hydrogen) are significantly 
softened by electron doping. Figure[5ja) and (b) show the 
softening of the hydrogen-related phonon modes, namely 
the A*, Bxui B^ g , and B 3g phonons. The softening of 
these hydrogen-related phonons is a result of the weak- 
ening of the hydrogen bridge bonds upon charge doping, 
which in turn should enhance hydrogen release. 

In-plane boron modes (A g and Bi g ) show more sub- 
tle doping dependence. The B\ g mode hardens signifi- 
cantly with electron doping at low doping levels, reflect- 
ing a strengthening of the boron network. When doping 
level exceeds 0.5 electrons per cell, the B ig mode soft- 
ens slightly. At this doping level, excess electrons start 
to occupy high energy antibonding states. Although the 
system is still dynamically stable, thcrmodynamically it 
may prefer to release hydrogen in exchange for stabiliz- 
ing the boron network. The Ai g mode, on the other 
hand, softens at low doping levels. Since this mode in- 
volves the stretching of both the B— B bonds and hydro- 
gen bridge bonds, the initial softening of this mode shows 
the combined effects of strengthening the boron network 
and weakening the hydrogen bridge bonds. Clearly, at 
small doping, the weakening of the bridge bonds domi- 
nate. 



Summary 

In summary, we predict a novel layered solid boron 
hydride structure (B2H2) and investigate its structural, 



5 




0.2 0.4 0.6 0.8 
Doping (electron/cell) 




0.2 0.4 0.6 0.8 
Doping (electron/cell) 



FIG. 5: Phonon energy as a function of doping levels 
for four hydrogen phonon modes at T for single-layer 
B2H2 a, the change in the phonon energy for the three 
higher energy H phonon modes BJg, Bi u , and A* g and 
b, for lower energy H phonon mode B3 g . 
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electronic and phonon properties using density functional 
based first-principles electronic structure methods. The 
structure consists of a graphene-like hexagonal boron net- 
work and bridge hydrogens which form 3c2e-like multi- 
center B— H bonds. The formation of multi-center B— H 
bonds is critical to stabilize the electron-deficient hexago- 
nal boron network. The phonon dispersion of single-layer 
B2H2 does not show any soft phonon modes, thus con- 
firming the dynamical stability of the proposed structure. 
We also investigate the effects of charging on the dynam- 
ical properties of B 2 H 2 . The softening of the hydrogen 
related phonon modes upon charging suggests a weaken- 
ing of the hydrogen bridge bonds. With the presence of 
extra electrons, the boron network becomes less electron- 
deficient therefore reducing the strength of the multi- 
center B— H bonds. These results suggest a charging- 
assisted hydrogen release mechanism. The boron net- 
work, however, remains stable upon charging. 
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